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PACS 91 . 30 . Px - Earthquakes 
PACS 89 . 75 . Kd - Patterns 

PACS 43 . 60 . Cg - Statistical properties of signals and noise 

Abstract. - We report remarkable similarities in the output signal of two distinct out-of- 
equilibrium physical systems - earthquakes and the intermittent acoustic noise emitted by crum- 
pled plastic sheets - Biaxially Oriented Polypropylene (BOPP) films. We show that both signals 
share several statistical properties including the distribution of energy, distribution of energy in- 
crements for distinct time scales, distribution of return intervals and correlations in the magnitude 
and sign of energy increments. This analogy is consistent with the concept of universality in com- 
plex systems and could provide some insight on the mechanisms behind the complex behavior of 
earthquakes. 



Understanding the underlying mechanisms that govern 
he complex spatio-temporal behavior of earthquakes is a 
stimulating challenge [TJ[2] . Concepts and methods from 
statistical physics has been largely applied to study earth- 
quakes, contributing to identify several patterns in seismic 
^) activity [3-11 . This approach also has been contributing 
to identify universal behavior in earthquakes - similarities 
between seismic records and the output signal of systems 
t-H in different research areas. 

' * For example, it has been reported that 7-ray events 
emitted by neutron stars and earthquakes share several 
• • distinctive statistical properties indicating tectonic activ- 
ity on neutron stars - 'starquakes' - analogous to earth- 
quakes on Earth [l2] . Another example is a reported anal- 
ogy between earthquakes and the Internet. Specifically, it 
has been found that two known empirical power laws for 
earthquakes - the Omori law and the Gutenberg-Richter 
law - hold also for the Internet (ping experiment). In this 
context, sudden drastic changes of the Internet time series 
are referred as 'internet quakes' 
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Earthquake pat- 
terns also can be observed in financial markets. It has been 
reported, for instance, that stock price fluctuations follow 
a power law distribution with exponent ~ 3 15 . This 



behavior is quantitatively similar to those found in earth- 
quakes (see refs. [8j[l3]), indicating an analogy between 



natural and financial earthquakes. For other examples of 
universal behavior in complex systems, see refs. 16 - 21 . 

Here, we compare earthquakes with the output signal 
of an out-of-equilibrium physical system - the intermittent 
acoustic noise emitted by crumpled plastic sheets. Some 
out-of-equilibrium physical systems emit crackling noises 
as a response to external conditions through events span- 
ning a broad range of sizes 22 - 29 . In particular, the sharp 



and intermittent noises emitted by some kinds of crumpled 
papers and similar materials - including plastic sheets - 
qualitatively remember earthquakes which arise when two 
tectonic plates rub each other. Starting from this qualita- 
tive picture, we search for a quantitative support for this 
analogy. Specifically, in a series of experiments we mea- 
sure the acoustic noise emitted by a crumpled plastic sheet 
- a Biaxially Oriented Polypropylene (BOPP) film - in re- 
laxation and compare these records with real data on the 
magnitude of earthquakes. We find that both processes 
exhibit several similar statistical properties. 

To quantitatively test this analogy, we consider real data 
on earthquakes obtained from the Northern California cat- 
alog for the period 1966-2006 30 . This seismic database 
contains ~ 435,000 records from one of the most active 
and studied geological faults on the Earth - the San An- 
dreas Fault. For each event, we calculate a measure of the 
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Fig. 1: Obtaining the energy intensity from the recorded acous- 
tic noise. Top: A sample of crumpled plastic sheet in relax- 
ation. Middle: The corresponding acoustic noise recorded (in 
arbitrary units). Bottom-left: Detail of a single event; and 
Bottom-right: the corresponding intensity of a single event, 
/ = cp 2 , where c is a constant and p is the sound pressure. 
The variable S, given by the area under the corresponding 
set of peaks, is related to the energy dissipated in the event. 
Because of the short time distance between events in the begin- 
ning of the relaxing process, it was not possible to discriminate 
individual events. Then, for each experimental dataset, we cut 
out the first 10 seconds in the recording. 



energy dissipated, E = exp(M), where M is the reported 
magnitude of the earthquake. 

We also perform a series of experiments in order to ob- 
tain measures of the energy dissipated by crumpled plastic 
sheets in relaxation process. First we crumpled a given 
sample of plastic sheet (of size 1.0m x 1.5m) into a com- 
pact ball. This procedure is similar to the common expe- 
rience of crumpling an unwanted sheet of paper prior to 
disposing of it [24]. In such conditions, the plastic sheet 
emit sound in discrete pulses of a variety of intensities for 
a relatively large time after released to relax (about 10 
minutes) . 

We record the sound emitted with a condenser micro- 
phone (Shure Microflex MX202W/N) positioned at 1 me- 
ter from the sample for five minutes and digitalized at fre- 
quency 8000 Hz. For the analysis, the noise was reduced 
by applying a cutoff filter for low frequencies. A single 
event is identified by a set of peaks with intensity larger 
than a threshold value. The start time of the single sound 
is taken as the time corresponding to the first sound with 
intensity larger than the threshold, whose value was cho- 
sen above the noise intensity. The end time corresponding 
to the time when the intensity becomes lower than the 
threshold for a time bigger than a given t c ~ 1 — 2ms (the 
characteristic length of single event). Figure 1 shows a 
typical example of the acoustic emission recorded and the 



corresponding energy intensity S obtained for a specific 
event. 

We first determine the probability distribution of the en- 
ergy intensity for earthquakes and crumpled plastic sheets. 
Figure 2a shows that both distributions are consistent 
with a power law decay, 



P(x) 



(1) 



with a ~ 3. This result, already known for earthquakes, 
suggests that both systems are self-organized into a scale- 
free state - there is not a typical scale for the dissipated 
energy. Observe that the power law exponent is quantita- 
tively similar for earthquakes and crumpled plastic sheets. 

To find information on the dynamics of energy dissipa- 
tion, we define energy increments as AE = E(i + r) — E(i) 
and AS = S(i + r) — S(i), where E(i) and S(i) are pro- 
portional to the energy of i-th event. The distributions of 
energy increments, for different values of the time scale r, 
are shown in Figures 2b (earthquakes) and 2c (crumpled 
plastic sheets). All curves are symmetrical, very peaked 
and have wings larger than expected for a normal process. 
In both cases, data for distinct time scales collapse onto a 
single curve indicating that the distribution of energy in- 
crements exhibits a common functional form for all time 
scales in the range considered. We also shuffled the origi- 
nal series and then calculated the distribution of AE and 
AS again, but no significant changes were observed. This 
result may indicate no correlations or weak correlations in 
the time organization of energy intensities. 

Assuming a given variable x following a power law dis- 
tribution with exponent a (see eq. 0) and no correlation 
between two events (a first approximation), the probabil- 
ity distribution for the increments Ax = x(i + r) — x(i) is 
given by P(Ax) = K dx dx' {xx') ot 5{x' — x — e) = 
K dx[x(x + |e|)] _Q; , where K is a normalization con- 
stant and e is a small positive value to avoid divergence in 
x = 0. The integration leads, for real and positive a, to 
the normalized probability density function (PDF) 
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where F-j 2 is the confluent hypergeometric function [§]. 
This PDF is shown in Figures 2b and 2c in comparison 
with real data. Notice that both curves are given by eq. 
Q with the same parameters. The good adjustment to 
the data indicates that the distribution of energy incre- 
ments exhibits a common shape for both systems. 

The non-Gaussian behavior of the distributions of en- 
ergy increments, shown in Figures 2b and 2c, also can 
be characterized by q-Gaussian distributions - typical in 

In fact, 



Tsallis statistics 



31-33 



it has been reported 
that eq. Q can be very well reproduced by means of 
q-Gaussians, whose values of q are related with the pa- 
rameter a 1 8 . In the range considered, a q-Gaussian dis- 
tribution, with q ~ 1.75, practically coincides with the 
curves shown in Figures 2b and 2c (solid lines). For q > 1, 
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Fig. 2: Analysis of the probability density function (PDF) of data, (a) Probability density, P(E /(Je) and P(S/as), of normalized 
energy for earthquakes (circles) and crumpled films (squares), a is the standard deviation calculated over all records in a given 
series. The curves are shown vertically shifted for clarity. The solid lines are power laws given by eq. ([!]), with exponent a = 3.2. 
(b) Probability density of normalized increments, P(AE/<jae), for time scales r = 1, 10, 100. The solid line corresponds to 
equation ^ with a = 3.2. (c) Probability density P(AS/gas), for time scales r — 1, 10, 100. As well as in b, the solid line is 
given by equation Q with a — 3.2. (d) Probability density of normalized return intervals, P(R//j.r), for earthquakes (circles) 
and crumpled films (squares). In both cases, R is calculated within subseries of size 1016. /hr is the average of R in a given 
subseries. The threshold is R c = 2. The solid line is an exponential distribution P(r) = exp[— r]. 



the tails of a q- Gaussian decreases as a power law with ex- 
ponent /3 = 2/(q — 1). This result indicates that the tails 
of the distribution of energy increments follows a power 
law behavior, 

P(Ax)~Ax-P, (3) 

with fj ~ 2.7 for both systems. 

Another way to characterize the dynamics of the output 
signal of a given system is to analyze the return interval 
series. The return intervals R are defined as the interval 
between events that exceed a certain threshold R c . We ob- 
tain R from the normalized energy series - with elements 
E/ge an d S/gs, where g is the standard deviation. Fig- 
ure 2d shows the distribution P(r) for earthquakes and 
crumpled plastic sheets, where r — Rj \ir and \ir is the 
average of R in a given subseries of size 1016 (the typical 
size of a given experimental record for a crumpled sheet). 
For comparison, we also show the exponential distribution 
P{r) = exp[— r]. Observe that the distribution of return 
intervals for earthquakes and crumpled films share a com- 



mon shape. We perform a parallel analysis for several val- 
ues of the threshold R c but no significative changes were 
observed. The exponential behavior found in the distri- 
bution of return intervals suggests no correlation or weak 
correlations in the energy series of both systems. 

Next, we investigate fractal properties in the output sig- 
nal of the systems. For a given time series x(i), the auto- 
correlation function is defined as C(r) = (x(i)x(i + r)) — 
(x(i) 2 ). In order to reduce fluctuations in C(r) it is com- 
mon to obtain the root mean square fluctuation function 
F(r), such that F(r) 2 = £U £J=i ^0" - 1) [38 . The 
net displacement after r steps is y(r) = Ym=i^v) an ^ 
the root mean square fluctuation is defined as F{r) = 
V (Ay(r) 2 } - (Ay(r)} 2 , where Ay(r) = y(r + r) - y(r ). 
For fractal series, F{r) follows a power law behavior, 
F(t) ~ r h , where h is the scaling exponent which quan- 
tifies the degree of correlations. When h > 0.5 (h < 0.5) 
the series is long-range correlated (anti-correlated). Un- 
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Fig. 3: Analysis of temporal correlations in the data. F av (n) is an average of the DFA fluctuation function F(n) calculated 
for several subseries. F^ l {ri) is the average DFA fluctuation function obtained from integrated subseries. We investigate F in 
the range 5 < n < 100. We also show F av and F™ 1 obtained from shuffled series (open circles). In all cases, the exponent 
h is obtained by least square linear fits to the data. As expected, h c± 0.5 for all shuffled series, (a) F av (n) calculated from 
normalized energy series for earthquakes, E/a(E). The slope gives h = 0.57. (b) F av (n) calculated from normalized energy 
series for crumpled plastic sheets, S/a(S). The linear fits gives h = 0.55. (c) F av (n) calculated from normalized magnitude series 
for earthquakes, \AE\/a(\AE\), giving h = 0.66. (d) F av (n) calculated from normalized magnitude series for crumpled plastic 
sheets, | AS|/cr(|AS|). In this case, h — 0.63. (e) F^in) calculated from integrated sign series for earthquakes, sign(Ai^). The 



slope is 1.30, giving h = 0.30. (f) 
slope is also 1.30, giving h = 0.30. 



calculated from integrated sign series for crumpled plastic sheets, sign(AS). The 



correlated series present h = 0.5. Short-range correlations 
also may exhibit h = 0.5. 

For nonstationary records, it is common to apply de- 
trended fluctuation analysis (DFA) [34j[35] to investigate 
correlations in the data. For fractal series, the DFA root 
mean square fluctuation, F(n), also follows a power law 
behavior, 

F(n) ~ n h , (4) 

where n is a time scale. Here we apply DFA method in 
order to quantify temporal correlations in the data. 

The typical size of a given time series in the experiment 
of crumpled plastic sheets is ~ 1016 (20 samples). In or- 
der to perform a parallel analysis, we partition the energy 
series for earthquakes in subseries of size 1016 (428 sam- 
ples). For each subseries we obtain the DFA root mean 
square fluctuation, F(n), and perform an average of F(n) 
over all subseries - obtaining F av (n). Because the size 
of a typical subseries, we investigate F(n) in the range 
5 < n < 100. Figures 3a and 3b shows F av (n) for the 
energy series for earthquakes and crumpled plastic sheets 
- E and S. We find h ~ 0.55 for both records suggest- 
ing weak correlations in the data. As expected, we find 
h ~ 0.5 for shuffled series. 

Starting from the sequence of energy increments, we 



also obtain two sub-series: magnitude of energy incre- 
ments - \AE\ and |A5| - and sign of energy increments 
- sign[A£?] and sign [AS]. The function sign [Ax] assumes 
the values —1, or 1 if the increment Ax is negative, null 
or positive, respectively. For details of the magnitude- 
sign decomposition approach, see refs. 36 37 . Figures 3c 
and 3d shows F av (n) for the magnitude series of energy 
increments for earthquakes and plastic sheets. For both 
records h c± 0.65 indicating long range correlations in the 
data. Figures 3e and 3f shows F^(n) - the average DFA 
fluctuation function obtained from integrated series - for 
the sign series of energy increments for earthquakes and 
plastic sheets. This previous integration is necessary in 
DFA method when h < 0.5. For both records we find 
h ~ 0.30 indicating anti-correlations in the data. As ex- 
pected, shuffled magnitude and sign series exhibits h ^ 0.5 
indicating uncorrelated behavior. Notice the quantitative 
agreement between the values of h for both systems. 

The analysis reported here indicates remarkable simi- 
larities between two distinct out-of-equilibrium physical 
systems, providing a quantitative support for the analogy 
between earthquakes and crumpled films. Specifically, we 
show that for both signals i) the distribution of energy 
follows a power law with exponent a ~ 3; ii) the dis- 
tribution of energy increments exhibits a common non- 
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Gaussian shape in the range 1 < r < 100, with power 
law tails with exponent j3 ~ 3; iii) the distribution of 
return intervals follows an exponential behavior; iv) the 
DFA power law exponent is h ~ 0.55 for energy series, 
h c± 0.65 for magnitude series of energy increments and 
h ~ 0.30 for sign series of energy increments in the range 
5 < n < 100. These findings are consistent with the hy- 
pothesis that earthquakes and crumpled plastic sheets may 
be driven by common underlying mechanisms. 

The nature of both processes analyzed here also presents 
analogies. It has been pointed that the energy stored in 
a crumpled material is originated in the buckling process 
while the film is crumpled, and it is mainly concentrated 
in the formed ridges 39,40 . The nonequilibrium behavior 
observed in the relaxation process can be understood as a 
consequence of the frustration in the crossed ridges, char- 
acterizing a stress [41]. Moreover, it is known that Earths 
crust can also exhibit buckling under viscous stresses on 
its layers [40] . 

Some phenomenological models, as the epidemic- 
aftershock sequence model (ETAS) 



43 , the 



type aftershock sequence model (ETAS) [42j 
continuous-time random walk models (CTRW) [43 and 
the Olami-Feder-Christensen model (OFC) 8 : 44 45 try to 
incorporate the main properties of the complex spatiotem- 
poral behavior of earthquakes. Since the experiments with 
crumpled plastic sheets are simple and reproducible, they 
may be used as an additional data source to compare with 
artificial data. We hope that this analogy could provide 
some insight on the mechanisms behind the complex spa- 
tial and temporal behavior of earthquakes. 
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